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ABSTRACT 

A standard perturbation expansion in the a tom- radiation field inter- 
action is used to calculate the two photon excitation cross section for 
the Is^ 2s^ 2p** ^P to Is^ 2s^ 2p^(**S) 3p ^P transition In atomic oxygen. 
The summation over bound and continuum intermediate states is handled by 
solving the equivalent inhomogeneous differential equation. We find 
that our exact summation results differ by a factor of 2 from a rough 
estimate obtained by limiting the intermediate state summation to one 
bound state. Higher order electron correlation effects are also examined. 


♦National Research Council Research Associate 


I. INTRODUCTION 

Atmospheric remote sensing systems are currently being developed 
which will take advantage of the remarkable properties of high-energy 
lasers.^ With detection and ranging systems based on optical frequencies 
one may obtain information on atmospheric pressure, temperature, and 
water vapor altitude profiles as well as particulate number densities. 

For certain projected applications^ the height distribution of trace 
atomic and molecular gases in the stratosphere may be obtained by knowing 
their various two-photon excitation cross sections. 

Recent multi photon work^ has shown that the validity of standard 
perturbation theory extends to quite strong photon intensities. For 
two photon absorption the sum over intermediate bound and continuum 
states found in the second order expression for the cross section may 
be handled by solving the equivalent inhomogeneous differential equation. ‘‘*5 
This technique has been applied to two photon ionization of hydrogen®*^ 
and the alkali atoms®"^°, as well as recent studies of double photo- 
ionization^^ and the hyperfine interaction. ^2 

In this paper we use perturbation theory to calculate the two 
photon excitation cross section for the Is^ 2s^ 2p** to Is^ Zs^ 2p^ (‘♦S) 

3p transition in atomic oxygen. We find that our exact sumnation 
results using the differential equation technique differ by a factor of 
2 from a rough estimate obtained by limiting the intermediate state 
summation to one bound state. Certain higher order electron correlation 
corrections to the second order expression are also investigated. In 
Sec. II we present the perturbation theory for the two photon excitation 
of an acorn; Sec. Ill contains the results of tf.e oxygen calculation; and 
Sec. IV is a brief summary. 
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II. THEORY 


The photoexcitation cross section for an atom may be calculated using 
t1n« dependent perturbation theory, The non-relativistic Hamiltonian for 
an N electron atari is given by 




= t Hi ) 
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The single particle potential Vj is generally the Hartree-Fock potential, 
while V.. is the Coulomb interaction between pairs of electrons. In Eq. 

(2) m is the electron mass, e is the magnitude of the electronic charge, 
and Z is the atomic number. 

The unperturbed Hamiltonian H^ for the atom-radiation field system is 
obtained by combining of Eq. (2) with the Hamiltonian for the quantized 
electromagnetic field. jhe eigenstates |*j> of Hq are products of an 
electronic state |j> and a radiation state with n^ photons per mode x such 
that n = En . The eigenenergy E*? equals E. + l n.fi<u. , where E. is the un- 

X^ j JX''X j 

perturbed energy of the jth electronic state and is the angular frequency 
of the xth mode. 

The perturbation Hamiltonian H^ is given by 

Hi = , (4) 
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and the annihilation and creation operators for the quantized vector 

1 3 

potential A(f) may be -dfined by 
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The solution v(t) of the time-dependent Schrodinger equation with 
Hamiltonian H = Hq + 1s given by^** 
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and the parameter y allows v(t) to be adiabatically developed from the 
Initial state |*q>. By projecting v{t) of Eq. (9) onto the final state 
> and then carrying out the time integrations in Eq. (10), the transition 
rate R is found to be^^ 


fv ^ I *' ' 

In a many-body diagrammatic expansion evaluation of Eq. (12), the summation 

L means only linked diagrams are to be included, 

For a single beam of light the two photon transition rate for 
excitation may be evaluated from Eq. (12) by considering a radiation 
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field containing only one mode.^^ In the dipole approximation the 
lowest order expression for R with m = n-2 is given by 


R = F"- I ^ S(a«,-Ay)) , 

where F is the photon flux (in photons/cm^ sec) at frequency ut Aw ® 
and I represents a sum over bound states and an integration over the 
continuum. If we assume that the final state has some finite lifetime 
T^, we may replace the delta function in Eq. (13) by^® 


S(au)-Ai>)) 


i 
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We may now define a generalized cross section 8 (in cm** sec) for two 
photon absorption as 

^ • 05) 

At the peak frequency of the Lorentz distribution, 55 = _jw, the generalized 

2 

cross section in the velocity form of the dipole operator is: 


= itnVK* 




(17) 


I < f| |,^?jln>lX^»l||,^V'Jo> I ( 16 ) 

^ E,-E.> + «i5 

while the length form yields: 

$L * l (»iT*e*S *’ Tt \ t 

c* I •» ^ 

The peak cross sections of Eqs, (16) and (17) will not depend on 
the actual laser bandwidth. A total two photon cross section, how- 
ever, must be integrated over the frequency distribution of the 
particular light source. We should also point out that our cross 
section equations have been derived assuimning that the laser pro- 
duces photons in a purely coherent state. Any photon statistics effects 
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or bunching of the photons will Increase the probability for absorption 
by as much as a factor of 2.^ 

Higher order electron correlation and radiation field corrections 
will add to the terms Inside the absolute squares of Eqs. (16) and (17). 
For exact atomic wave functions the velocity and length forms of the 
dipole operator are of course equivalent, but since we use a perturbation 
expansion for they may differ by significant amounts in lowest 
order. 
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III. CALCULATION AND RESULTS 

In this section we epply the perturbation theory of Sec. II to the 
calculation of the two photon excitation cross section for the Is^ 2$^ 

2p** 3p to Is^ 2s^ 2p3 (*♦$) 3p 3p transition In atomic oxygen. The experi- 
mental term separation (weighted average over J levels) Is 10.98 eV^^, 
thus a laser energy of 5.49 eV (or 2258 %) is needed for two photon excita- 
tion. Unless otherwise indicated we will assume that the Incident radiation 
is linearly polarized. 

To obtain the lowest order generalized cross section we must 
evaluate Eqs. (16) and (17). In the case of the 2p** ground state of 
oxygen, we calculate 6 for * 0 and also for Mj_ * + 1 and then average 
over these v lues. The single particle core orbitals for the ground 
state were calculated numerically in the Hartree-Fock approximation . 20 
In lowest order the intermediate virtual states allowed by the dipole 
selection rules are Is^ 2s^ Zp^ (*♦$) ns Is^ 2s^ 2p^ (“^S) nd 

and Is^ 2s 2p‘‘ (*♦?) 3p ^0°, ^P®, We used the Silverstone^^-Huzinaga^^ 

frozen core potential to construct the orthogonal excited state spectrum. 

This potential has the property of preserving the Hartree-Fock operator 

N-l 

for the core orbitals while at the same time providing a realistic V 
type operator for the excited states. We assume all single particle states 
have the product form P„£(r) where Y^(0»4) is a 

spherical harmonic andXj 1s the spin function. The potential for 
the radial orbitals (both bound and continuum) of the 2p^ (‘*5) ns 

intermediate state is then given by: 

Vs«(^o) - - kA 

t i K‘p , 
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where 
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and r^« tnin (r^, r) and * max Cr^, r). The potential for the radial 
P^d(>^) orbitals of the 2p^ (^S) nd 3D® Intermediate state Is given by 

= «iJiS ” ? Kis + ^Jjls ' i { 72 ) 

^ 3 Jip r ^ Kip f ^ Kap • 


As a first approximation the P^pCr) orbital of the 2p3 (‘♦S) 3p final 
state of oxygen, with single particle energy E 3 p, is obtained by solving 


( Ti - ^ 0 , (23) 
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For the cross section b(+ 1) the dipole selection rules for determinantal 
states allow only contributions from the 2p^ (**S) nd intermediate states 
in lowest order. After performing the angular integrations, we obtain the 
following radial form of Eq, (16) for the 2p** ^P to 2p^ (‘♦S) 3p 
transition in oxygen: 

' l UirVl * T. r fT 

c* L 
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where 


J 


(27) 




1 Pn'VlO Dv PnlU) 


^ ^n£ II Oy 11 4p^ (28) 

nl £ap- £m + "tliO 

and Oy * d/dr + ~ One should note that 

the Intermediate sum Is now present In the expression for the perturbed 
bound state function (^^(n-f.r) of Eq. (28). For the cross section b( 0) 
the dipole selection rules allow contributions frcxn all three types of Inter- 
mediate states. The radial form of Eq. (16) Is thus; 


^(#1 = t.r ♦ f? < 3 pii 0, 

L'l" ^ (29) 

4 [T < isilP.IUpXipllD.lia» > 1 . 

V343 ^ ^ 1\u5 , 

The radial matrix elements for Bj_(i.l) and $j^(0) derived from Eq. (17) 
are the same as those found In Eqs. (26) and (29) except that Is 
replaced by ^ r and (rvf,r) is replaced by 4 >j^(n£,r). 

The perturbiid bound state function <p^(nl,r) of Eq. (28) may be 
found by solving the Inhomogeneous differential equation***^: 


( \ - aV;,(i) 4 if.) '^,("t,r) ^ Dv Riplf) , (30) 

where C 2 p + w. The exper1n«ntal energy of C 2 p = - 13.62 eV^^ is 
used for functions (jiy(n^,r) which represent 2p5(**S) nt ^L® virtual 
excitations. Similar equations for perturbed functions ♦j_(n£,r) were 
also solved. The accuracy of our solutions was checked by computing 
cross sections for the two photon ionization of argon using perturbed 
bound state functions and comparing them with results obtained using the 
radial Integration truncation method. two methods agree to three 
significant figures In both the length and velocity forms at energies 
away from resonances. 
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The lifetime of of the 2p^ (‘♦S) 3p ^p final state of oxygen may 
be taken as the natural decay time for single photon emission to the 
2p3 ('♦S) 3s state. This will only be true when the laser pulse time 
Is such that T|^ > For laser Intensities currently envisioned 

In atmospheric remote sensing work? the photoionization rate from the 2p^ 

(•♦S) 3p 3p final state will be much smaller than the natural decay rate. 

The handbook^^ value of * 3.57 • 10’® sec (accuracy to 10%) is used 
for all subsequent calculations. If we use Eqs. (A2) and (A3) of the 
APPENDIX along with lowest order that * 3.31 • 10"®sec 

In the velocity form and ** 2.97 * 10*® sec In the length form. 

The lowest order Hartree Fock (HF) results for the generalized cross 
section 0 for the 2p** ®P to 2p® (‘♦S) 3p ®P transition In oxygen are given 
In both the length and velocity forms in Table I. Due In part to the 
Inclusion of contributions from the 2p® (**S) ns Intermediate states, 
b( 0) is seen to be a factor of 20 to 30 larger than b(+ 1). For b( 0) the 
velocity form is 1.4 times larger than the length form, while for 6(+l) 
the two forms almost agree. If we transform our wave function representa- 
tion according to: 

Y(lsjm,^ Z « (31) 

rt^Pis 

where {/l'^2^3 is the usual Cl ebsch-Gordan coefficient, we may construct 
cross section expressions between fine structure levels using the results 
already obtained for b( 0) and ft(+l). The experimental fine structure wave- 
lengths^® and our lowest order cross section results are found In Table II. 

It should be noted that for two photon absorption we may have aJ » + 2, ♦ 1 
and 0. 
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The Importance of the complete sum over intermediate states In the 
results for the cross section a may be shown by comparing the matrix 

element ■ <3p|| 0y|| ay(ns)> with 

» <3p MlXIUs><istllXlli p> ^ (32) 

€4^ - €|s + IUa) 

where p35(^^) Is the radial orbital for the 20^ (‘♦S) 3s state. We 
found that (In atomic units) * 0.671 and « 0.443, while ® 12.8 
and * 19.2. Thus If we use only one intermediate state In the evalu- 
ation of 0(0). we would find a factor of 2 reduction In the velocity 
form while the length form would show a factor of 2 Increase. 

If we assunw Incident radiation that is circularly polarized, one 
must perform the angular Integrations In Eqs. (16) and (17) using 
different dipole selection rules. The lowest order Hartree Fock results 
for the 2p** to 2p3 (•♦$) 3p ^P transition In oxygen are then < * 

1.704 • 10“**^ cm** sec and <8j^>- 0.935 • 10"**^ cm** sec. 

Electron correlation effects may be Important In certain multiphotor 
processes In complex atoms. A many-body dlagramnatlc expansion of 
Eg. (12) was developed previously to handle electron correlations In 
the two photon Ionization of argon. x^e configuration Interaction R matrix 
method has been applied to the calculation of the third-order optical harmonic 
coefficient for helium. 2** We may also use the multi -configuration Hartree- 
Fock (MCHF) methodic to calculate electron correlation corrections for 
multi photon processes In a manner similar to that used In recent photo- 
ionization work. 25 

We therefore examine first the third order electron correlation 
corrections to the second order two photon cross section for the 2p** ^P 
to 2p^(**S) 3p ^P transition In oxygen. The complexity of the ground state 
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correlation problem in oxygen may be appreciated by studying Table IV 
of a correlation energy paper by Marchetti, Krauss and Weiss. 

The largest contribution comes from pair excitations of the 2p shell 
(43%). Since 2p^ may couple to give ^D, and an^ 3d^ may couple 
to g* i^e have quite a number of allowed ways 

to attain an overall symmetry. To get a feel for the strength of 
the ground state correlation correction, we chose the largest angular 
components among those 2p shell pair excitations given by Marchetti et 
al. The mixing coefficients of a fixed core MCHF calculation for the 
configurations 2 p 2 (352 15 ) 3 p^ 3P (Sd^ iQ) ^P and 2 p^ 3 p 

(3s3d^0) 3P are given in Table IIIA. The coefficients a^ of the three 
excited configurations are all quite small. The contribution of the 2p^ 
3p ( 3 d^ iD) configuration to the cross section 3 ^ (+ 1) of Eq. (26) 
is given explicitly by: 



where = t P-^^>~ K»<U| 3 A>< 3 p|| P« || 3 jl > . (34) 

*>4 £ip-e«i+'K'rf 

We may again solve for the perturbed bound state function (Sd^) using 
the inhomogeneous differential equation technique. This two step method 
of accounting for ground state pair excitations should be equivalent to 
solving the inhomogeneous pair equation used by Garpman et al .^2 Expressions 
for 0 y(O), 0 l(± 1) and 3 l(^) have forms similar to Eq. (33), 

The largest contribution to the electron correlation corrections 
for the 2p3 (*♦$) 3p ^p final state of oxygen should come from the mixing 
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of the 2p3 Intermediate couplings. This particular open shell effect 
was found to be quite strong in recent photo ionization work on chlorine. ^7 . 28 
The mixing coefficients of a fixed core MCHF calculation for the configura- 
tions 2p3 (‘*S) 3p 3P, 2p3 (^D) 3p 3P, and 2p3 (2p) 3 p 3p given in 
Table IIIB. The MCHF radial 3p function was found to differ only slightly 
from the 3p function obtained using Eq. (23). The contribution of the 2p^ 
(^D) 3p 3p configuration to the cross section {+ 1) of Eq. (26) is 
given explicitly by: 

= IkjtVif tr f IT 

L'/15 (35) 

- /r ^ M , 

VbO ViSoo J 

where the perturbed bound state functions (ns) and (nd) are obtained 
by using Eq. (30) with ~ - 16.94 eV.^^ Expressions for B^CO), B|_(+ 1) 
and B^ (0) also have forms similar to Eq. (35). 

The multi -configuration Hartree-Fock (MCHF) results for the generalized 
cross section B for the 2p** ^p to 2p^ (‘♦S) 3p ^P transition in oxygen are 
given in both the length and velocity forms in Table I. Cross sections 
between fine structure levels are given in Table II. We have included only 
the excited configurations 2p^ 3p (Sd^ iD) ^P, 2p3 ('’-d) 3p ^P and 2p^ (^P) 

3p 3P in our final results. The MCHF values for b( 0)* B(+ 1) and < 6 > 
are all lower than the corresponding Hartree-Fock values. We must note 
that the agreement between length and velocity has not improved. This 
indicates we have only touched the surface of a very complex electron 
correlation problem and the continuing difference between the length and 
velocity forms seems to be a good indicator o^ this fact. What we have 
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learned is that the result of a more complete solution of the problem 
will probably not differ from our lowest order Hartree-Fock values by more 
than a factor of 2, 

As one goes to stronger and stronger photon intensitiesB thn 
standard perturbation theory begins to break down. In this case the 
electric component of the radiation field is becoming comparable to the 
Coulomb field which binds the atom. From previous work^^ we estimate 
that radiation field corrections, such as the fourth order process 
involving the absorption of 3 photons and the emission of one, will 
begin to significantly affect our lowest order results at a flux F of 
10^^ photons/cm^ sec. At higher intensities one must use an exact 
summation of a whole series of higher order radiative terms^® or 
go to an alternative formulation of the problem,^ We should also note 
that long before this regime is entered (around 10^5 photons/cm^ sec), 
the Lorentz distribution of Eq. (14) will be modified by either a decrease 
in laser pulse time or photoionization from the excited state. 

This in turn will affect the excitation cross sections e of Eqs. (16) 
and (17). 
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IV. SUMMARY 

In this peiper we determined a value for a two (>hoton excitation cross 
section In oxygen of sufficient accuracy for laser atmospheric remote 
sensing applications. We used the Inhomogeneous differential equation 
technique to handle the sums over bound and continuum Intermediate states 
found in a standard perturbation theory development. Although the electron 
correlation corrections are found to be quite complex, their overall 
effect should not significantly alter the lowest order Hartree-Fock 
results. In the future we hope to extend our work to the study of 
certain two photon transitions In nitrogen, chlorine, and argon. 
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APPENDIX 


The various single photon emission and absorption transition rates 
may be obtained from Eq. (12) with m = n + 1. In the dipole approxima- 
tion the lowest order expression for the absorption rate is given by: 

I Fx , (ft, 

where Au = (E^ - E^l/fi. The induced emission rate has a form similar 
to Eq. (Al) except that Au = (E^ - E^)/fi. For spontaneous ©nission the 
flux F equals c/V, where V is the volume of the radiation cavity. If 
the summation j is replaced by an integration over u with the appropriate 
density of modes,^® we obtain for the spontaneous emission rates at 
uT = Aw: 




(A2) 

and 

i<4l2 . 

(A3) 


The absorption cross section a (in cm^) may be obtained from Eq. (Al) by 
considering an external radiation field containing only one irode and replacing 
the delta function as in Eq. (14). The cross section o at uli ~ A(i) i s then 



If we take for the natural lifetime for decay to |o>, we recover the 
well-known result^ ^ 



where is the statistical weight of state |1>. 

For two photon enlssion processes we find an induced emission that 
is proportional to F^, a spontaneous emission that is independent of 
photon flux, and cross terms that are linear in F. The spontaneous 
emission rate may also be obtained from Eq. (12) by replacing the 
sumnatlons over modes and X 2 by integrations over frequencies, giving: 

Jo 


E.-t. 


^ (A7) 


M V 


> 75 - 
-<U) 



where u = - oj. In thermal equilibrium the level populations of states 

|o> and |f> are related by Boltzmann's law. Equating two photon absorption 
and emission processes at each frequency w, we may again recover Planck's 
law for the frequency distribution of the photon flux. 
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TABLE I: The Two Photon Cross Section b for the 2p** to 

2p3 {‘♦S) 3p Transition in Oxygen (in units of 10"**^ cm** sec) 


Approximation 

b(o) 

6(+l) 

<B > 

HF velocity 

6.538 

0.232 

2.334 

HF length 

4.729 

0.234 

1.732 

MCHF velocity 

6.330 

0.155 

2.213 

MCHF length 

4.566 

0.174 

1.638 


* < 8 > - [Z6(il) + 8(0) ]/3 
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TABLE II: The Fine Structure Components of the Two Photon Cross Section 
a for the 2p** to 2p^ (‘♦S) 3p ^P Transition In Oxygen (In 
units of cm** sec) 


^1 

'^f 

x(X) 



.MCHF 

-MCHF 

2 

0 

2256.5390 

0.192 

0.127 

0.200 

0.132 

2 

2 

2256.5430 

1.711 

1.319 

1.562 

1.211 

2 

1 

2256.5574 

0.431 

0.286 

0.451 

0.296 

1 

0 

2260.5758 

0.000 

0.000 

0.000 

0.000 

1 

2 

2260.5798 

0.718 

0.477 

0.751 

0.493 

1 

1 

2260.5940 

1.616 

1.256 

1.462 

1.146 

0 

0 

2262.3328 

1.377 

1.097 

1.21V 

0.980 

0 

2 

2262.3368 

0.957 

0.636 

1.001 

0.657 

0 

1 

2262.3510 

0.000 

0.000 

0.000 

0.000 
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TABLE III: Multi -Configuration Hartree-Fock Calculations for Oxygen 


Ground State 

Configuration 

Mixing Coefficients 


= 0.9989 

2p2 3p (3 s 2 IS) 3p 

a^ - 0.0191 

Zp 2 3p (3d2 ID) 3p 

a^ * 0.0398 

Zp2 3p (3s 3d ID) 3p 

84 * 0.0143 

Final State 

Configuration 

Mixing Coefficients 

2p3 (‘»S) 3p 3p 

bj = 0.9958 

2p3 (2 d) 3p 3p 

bg = 0.0862 

2p3 (2P) 3p 3p 

b 3 » - 0.0323 
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